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Abstract 

Stochastic network calculus is a newly developed theory for stochastic ser- 
vice guarantee analysis of computer networks. In the current stochastic net- 
work calculus literature, its fundamental models are based on the cumulative 
amount of traffic or cumulative amount of service. However, there are net- 
work scenarios where direct application of such models is difficult. This paper 
presents a temporal approach to stochastic network calculus. The key idea is 
to develop models and derive results from the time perspective. Particularly, 
we define traffic models and service models based on the cumulative packet 
inter-arrival time and the cumulative packet service time, respectively. Rela- 
tions among these models as well as with the existing models in the literature 
are established. In addition, we prove the basic properties of the proposed 
models, such as delay bound and backlog bound, output characterization, 
concatenation property and superposition property. These results form a 
temporal stochastic network calculus and compliment the existing results. 

Keywords: Stochastic network calculus, max-plus algebra, min-plus 
algebra, stochastic arrival curve, stochastic service curve, performance 
guarantee analysis, delay bound, backlog bound 
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1. Introduction 



Stochastic network calculus is a theory dealing with queueing systems 
found in computer networks [8] 15] 20] 22] . It is particularly useful for ana- 
lyzing networks where service guarantees are provided stochastically. Such 
networks include wireless networks, multi-access networks and multimedia 
networks where applications can tolerate a certain level of violation of the 
desired performance 
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Stochastic network calculus is based on properly defined traffic models 



20] 22] 27] 34] 37] and service models 20] 22] . In literature, it is typical to 



model the arrival process by a stochastic arrival curve and the service pro- 
cess by a stochastic service curve. The arrival curve provides probabilistic 
upper bounds on the cumulative amount of arrival traffic whereas the service 
curve lower bounds the cumulative amount of service. In this paper, we call 
such models space- domain models, from which extensive results have been 

(P.l) Service 
Output Charac- 



[22 



derived. There are five most fundamental properties 
Guarantees including delay bound and backlog bound; (P.2 
terization; (P.3) Concatenation Property; (P.4) Leftover Service; (P.5) Super- 
position Property. Examples demonstra ting the necessity and applications of 
these basic properties can be found 

However, there are many open challenges for stochastic network calculus, 
making its wide application difficult 22J. One is to analyze networks where 



users are served probabilistically. For example, in wireless networks, a wire- 
less link is error-prone and consequently retransmission is often adopted to 
ensure reliability. In random multi-access networks, random backoff and re- 
transmission are used to deal with contention and collision. To apply stochas- 
tic network calculus to analyze such networks, it is fundamental to find the 
stochastic characterization of the service time that provides successful trans- 
missions for the user. However, direct application of existing space-domain 
models, which are built on the amount of cumulative service, is difficult. 

This paper aims to rethink stochastic network calculus to address some 
of the challenges in current stochastic network calculus literature, such as 
the analysis of error-prone wireless channels and/or contention-based multi- 
access. To be specific, we present a temporal approach to stochastic network 
calculus. The key idea is to develop models and derive results from the 
time perspective. We define traffic models and service models based on the 
cumulative packet inter-arrival time and the cumulative packet service time, 
respectively. In this paper, we shall call such models time-domain models. 
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In addition to their easy use in network scenarios discussed above, the basic 
properties are also investigated based on the proposed time-domain models. 
Moreover, relations among the proposed time-domain models as well as with 
the corresponding space- domain models are established, which provide a tight 
link between the proposed temporal stochastic network calculus approach and 
the existing space domain stochastic network calculus approach. This gives 
increased flexibility in applying stochastic network calculus in challenging 
network scenarios. 

The structure of the paper is as follows. In Section [2] we first introduce 
the notation and the system specification, followed by the review of the rel- 
evant results of stochastic network calculus. Section [3] defines the network 
calculus models in the time-domain and explores the model transformations. 
Four fundamental properties are thoroughly investigated in Section HJ The 
relevant discussion reveals the reasons of establishing the model transforma- 
tions in Section [3j In Section |5l we conclude the paper and discuss the open 
issue. 

2. Network Model and Related Work 

This section specifies the network system and reviews mathematical pre- 
liminaries for the analysis in the following sections. A brief overview on 
stochastic network calculus of particular relevance to this paper is presented 
as well. 

In this paper, we make the following assumptions unless stated otherwise. 

• All packets have the same length. 

• A packet is considered to be received by a network element when and 
only when its last bit has arrived to the network element. 

• A packet can be served only when its last bit has arrived. 

• A packet is considered out of a network element when and only when 
its last bit has been transmitted by the network element. 

• Packets arriving to a network element are queued in the buffer and 
served in the FIFO order. All queues are empty at time 0. 

• All network elements provide sufficient buffer space to store all incom- 
ing traffic and are lossless. 
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2.1. Notations and System Specification 

We use P(n), a(n), d(n) and 5 n , to denote the (n + l)th packet entering 
the system, its arrival time to the system, its departure time from the system 
and its service time provided by the system, respectively, where n — 0, 1, 2, .... 

• From the temporal perspective, an arrival process counts the cumula- 
tive inter-arrival time between two arbitrary packets and is denoted by 
T(m, n) = a{n) — a(m) for any < m < n. Note T(n, n) = 0. 

• A service process describes the cumulative service time received be- 
tween two arbitrary packets and is denoted by A(m,n) = Ylk=m^k f° r 
any < m < n. Note A(n, n) = 5 n . 

In the time-domain, the system backlog and system delay are defined 
below, respectively. 

Definition 1. The system backlog at time t > is denoted by B(t): 

B(t) < inf {/ > 0, sup{n > : a(n) < t} : d(n -l)< a{n) } . (1) 

The delay that packet P{n) experiences in the system is denoted by D{n): 

D{n) = d{n)-a{n). (2) 

Moreover, the time that packet P(n) waits in queue is denoted by W(n): 

W(n) = D(n) - 5 n . (3) 

The following function sets are often used in this paper. 
The set of non-negative wide-sense increasing functions is denoted by J 7 , 
where for each function /(•), 

F={f(-):V0<x<y,0<f(x)<f(y)} 

and for any function /(•) G J 7 , we set f(x) = for all x < 0. 

We denote by T the set of non-negative wide-sense decreasing functions 
where for each function /(•), 

?={f(-):W<x<y,0<f(y)<f(x)} 

and for any function /(•) G we set f(x) = 1 for all x < 0. 
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We denote by Q a subset of T, where for each function /(•) G Q, its 
nth-fold integration, denoted by f( n '(x) = (f°° dy) n f(y), is bounded for any 
x > and still belongs to Q for any n > 0, i.e., 



For ease of exposition, we adopt the following notations in this paper: 



In addition, the ceiling and floor functions are used in this paper as well. 

• The ceiling function \x~\ returns the smallest integer not less than x. 

• The floor function \_x\ returns the larget integer not greater than x. 

2.2. Mathematical Basis 

An essential idea of (stochastic) network calculus is to use alternate al- 
gebras, particularly min-plus algebra and max-plus algebra jij], to transform 
complex non-linear network systems into analytically tractable linear sys- 
tems. To the best of our knowledge, the existing models and results of 
stochastic network calculus mainly based on min-plus algebra that has basis 
operations suitable for characterizing the amount of cumulative traffic and 
service. As a result, these models focus on describing network behavior from 
the spatial perspective. Max-plus algebra is suitable for arithmetic opera- 
tions with cumulative inter-arrival times and service times. Consequently, 
network modeling from the temporal perspective more relies on max-plus 
algebra. In the following, we review the basics of both min-plus algebra and 
max-plus algebra. 

In min-plus algebra, the 'addition' operation represents infimum or min- 
imum when it exists, and the 'multiplication' operation is +. The min-plus 
convolution of functions /, g G J 7 , denoted by <8>, is defined as 



where, when it applies, 'infimum' should be interpreted as 'minimum'. The 
min-plus deconvolution of functions /, g G J 7 , denoted by 0, is defined as 




+ — 



max[x,0] and [x]x = min[x, 1]. 



(f®g)(t)=m£{f( 8 )+g(t- 8 )} 

0<s<t 



(f0g)(t)=sup{f(s + t)-g(s)} 
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where, when it applies, 'supremum' should be interpreted as 'maximum'. 

In the max-plus algebra, the 'addition' operation represents supremum 
or maximum when it exists, and the 'multiplication' operation is +. The 
max-plus convolution of functions /, g G J 7 , denoted by £g>, is defined as 

(f®g)(n)= sup {/(m) + g{n-m)} 

0<m<n 

where, when it applies, 'supremum' should be interpreted as 'maximum'. 
The max-plus deconvolution of functions f,g G J 7 , denoted by 0, is defined 

as 

(f@9)(n) = inf {f(n + m) - g(m)} 

where, when it applies, 'supremum' should be interpreted as 'maximum'. 
The max-plus convolution is associative and commutative [5|]. 

• Associativity: for any gi,g 2 ,g 3 e J 7 , (#i<g># 2 )<i>#3 = 9\®{g2®gz)- 

• Commutativity: for any <7i,<?2 G J 7 , gi®g 2 = g-i®9\- 

2.3. State of The Art in Stochastic Network Calculus 

The available literature on stochastic network calculus mainly focuses 
on modeling network behavior and a naly zing network performance from the 
spatial perspective QQHHHHHHHH- We call the corre- 
sponding models and results space-domain models and results in this paper. 

In order to characterize the arrival process of a flow from the spatial 
perspective, let us consider the amount of traffic generated by this flow in 
a time interval (s,t], denoted by A(s,t). In the context of stochastic net- 
work calculus, the arrival curve model is defined based on a stochastic upper 
bound on the cumulative amount of the arrival traffic. Here, we only review 
one relevant space-domain arrival curve model, virtual- backlog- centric (v.b.c) 
stochastic arrival curve (SAC) [22]. 

Definition 2. (v.b.c Stochastic Arrival Curve) 

A flow is said to have a virtual-backlog-centric (v.b.c) stochastic arrival 
curve a(t) G T with bounding junction f{x) G T , if for all < s < t and all 
x > 0, there holds 



sup L4(s, t) — a(t 

0<s<t 



> x 



} < /(*)■ 



(4) 
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In stochastic network calculus, the service curve model is defined as a 
stochastic lower bound on the cumulative amount of service provided by the 
system. Two space-domain service curve models [22] are reviewed here. 

Definition 3. (Weak Stochastic Service Curve) 

A network system is said to provide a weak stochastic service curve (3{t) G 
T with bounding function g(x) G T for the arrival process A(t), if for all 
t > and all x > 0, there holds 

P{A®f3(t)-A*(t)>x}<g(x) J (5) 

where A*(t) denotes the cumulative amount of the departure traffic. 

Unlike the arrival curve, it is difficult to identify the service curve from 
(j3J) because it couples the arrival process, the service curve and the departure 
process. Thus we need a more explicit model to directly reveal the relation 
between the service process and its service curve such as the following model 

5j. 

Definition 4. (Stochastic Strict Service Curve) 

A network system is said to provide stochastic strict service curve (3(t) G 
T with bounding function g(x) G T , if during any period (s,t), the amount 
of service S(s,t) provided by this system satisfies, for any x > 0, 

P{S(s,t) < P(t- s) -x) < g{x). (6) 

Definition |4] is applied to 'any period' which implies both worst-case sce- 
nario and other scenarios. If we could determine a function 0(t) which makes 
Eq.Q hold under the worst-case scenario, then Eq.Q automatically holds 
under other scenarios as well. 

Based on the arrival curve and service curve models, five fundamental 
properties have been proved to facilitate tractable analysis. For example, 
they can be used to derive service guarantees including delay bound and 
backlog bound, characterize the behavior of traffic departing from a server, 
describe the service provided along a multi-node path, determine the arrival 
curve for the aggregate flow, and compute the service provided to each con- 
stituent flow. 
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• P.l: Service Guarantees (single-node) 

Under the condition that the traffic arrival process has an arrival curve 
a(t) with bounding function f(x) and the network node provides ser- 
vice with a service curve model /3(t) with bounding function g(x), the 
stochastic delay bound and stochastic backlog bound can be derived. 
Particularly, the backlog bound is related to the maximal vertical dis- 
tance between a(t) and j3(t); the delay bound is relevant to the maximal 
horizontal distance between a(t) and (3(t). 

• P.2: Output Characterization 

To analyze the end-to-end performance of a multi-hop path, one option 
is the node-by-node analysis approach. This approach requires being 
able to characterize the traffic behavior after the traffic has been served 
and leaves the previous node. The output process of a flow from a node 
can also be characterized by an arrival curve which is determined by 
both the arrival curve of the arrival process and the service curve of 
the service process. 

• P.3: Concatenation Property (multi-node) 

Network calculus possesses an unique property, concatenation property, 
which is also used to analyze the end-to-end performance but improves 
the results obtained from the node-by-node analysis. The essence of 
the concatenation property is to represent a series of nodes in tandem 
as a 'black box' which can be treated as a single node. The service 
curve of this equivalent system is determined by the service curve of all 
individual nodes along this path. 

• P.4: Superposition Property (aggregate flow) 

Flow aggregation is very common in packet-switched networks. If mul- 
tiple flows are aggregated into a single flow under the FIFO order, the 
aggregate flow also has an arrival curve which is the summation of the 
arrival curve of all constituent flows. 

• P.5: Leftover Service Characterization (per-flow) 

The leftover service characterization makes per-flow performance anal- 
ysis feasible under FIFO aggregate scheduling. The crucial concept is to 
represent all other constituent flows as an 'aggregate cross flow' which 
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can be characterized using an arrival curve. Then the service provided 
to the constituent flow of interest can also be described by a service 
curve which is determined by the service curve provided to all arrival 
flows and the arrival curve of the 'aggregate cross flow'. 



The superposition property of the v.b.c SAC [22] is reviewed here because 
it is relevant to the model transformation in the following content. 

Theorem 1. Consider N flows with arrival processes Ai, i = 1,...,N, re- 
spectively. If each arrival process has a v.b.c SAC ai(t) G T with bound- 
ing function f\(x) G T , then the aggregate arrival process has a v.b.c SAC 
a(t) G J 7 with bounding function f(x) G T , where 

N 

a(t) = ^ a i( f ) and f( x ) = h ® • ■ - ® Jn{x). 

i=l 

3. Time-domain Modeling and Transformations 

This section defines traffic and service models in the time-domain. Par- 
ticularly, traffic models are defined based on probabilistic lower bounds on 
the cumulative inter-arrival time between two arbitrary packets. Service 
models are defined in terms of the virtual time function and probabilistic 
upper bounds on the cumulative service time between two arbitrary packets. 
Moreover, we establish the transformations among these models as well as 
the transformation between the time-domain model and the space-domain 
model. 

3.1. Time-domain Traffic Models 

Consider an arrival process that specifies packets arriving to a network 
system at time a(n), n = 0,1,2,.... In order to stochastically guarantee 
a certain level of QoS to this arrival process, this arrival process should 
be constrained. By characterizing the constrained arrival traffic from the 
temporal perspective, we define an inter-arrival-time (i.a.t) stochastic arrival 
curve model. 

Definition 5. (i.a.t Stochastic Arrival Curve) 

A flow is said to have an inter-arrival-time (i.a.t) stochastic arrival curve 
A(n) G J 7 with bounding function h(x) G IF, if for any m, n > and x > 0, 
there holds 

P|o(m + n) - a(n) < [\(n) - x] + | < h(x). (7) 
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Eq.flTJ) indicates that function A(n) is a probabilistic lower bound on the 
cumulative inter-arrival time. The violation probability that the cumulative 
inter-arrival time is smaller than A(n) is bounded above by function h(x). If 
h(x) = for all x > 0, Eq.([7]) represents a time-domain deterministic arrival 
curve [lo[ which is a special case of the i.a.t SAC. 

Queueing theory typically characterizes the arrival process using the prob- 
ability distribution of the inter-arrival time between two consecutive cus- 
tomers: 



Comparing F(x) with Eq.([7|), we notice that Eq.((7]) gives a more general 
probability expression of the inter-arrival time between two arbitrary packets. 
From this viewpoint, F(x) is a special case of Eq.([7]). 

Example 1. 

Consider a flow of packets with fixed packet size. Suppose that packet 
inter-arrival times follow an exponential distribution with mean Then, 
the packet arrival time has an Erlang distribution with parameter (n, /x) 
where n denotes the number of arrival packets. For any two packets P(m) 
and P(m + n), their inter-arrival time satisfies, for x > 0, 



where y = - — x. Thus, the flow has an i.a.t SAC A(n) = -. 

The i.a.t SAC is simple but has limited applications. For example, con- 
sider a virtual single server queue (SSQ) fed with the arrival traffic which 
has an i.a.t SAC X(n) with bounding function h(x). Suppose that the virtual 
SSQ provides a constant service time A(l) for each packet. From Eq.Q, the 
waiting delay of P{n) experienced in the virtual SSQ is 



P{a{n) - a(n - 1) < x} = F(x). 





n-l 



k=0 



W(n) 



d(n) — a(n) — A(l) 



sup 





0<m<n 





10 



where a(m) is the beginning of the backlogged period within which packet 
P(n) is transmitted. Eq. (15|) is derived from the departure time given in 
Eq.f JTB"]) . Eq.(j9]) is called the virtual- waiting- delay property. It is difficult 
to compute the virtual- waiting-delay from Eq.([7]). When investigating the 
performance guarantees such as delay bound and backlog bound in Section 
14.11 we f &ce the similar difficulty. 

In order to deal with the difficulty of computing the virtual-waiting-delay, 
we define another stochastic arrival curve model based on Eq.(j9]). 

Definition 6. (v.w.d Stochastic Arrival Curve) 

A flow is said to have a virtual- waiting- delay (v.w.d) stochastic arrival 
curve X(n) G J 7 with bounding function h(x) E T ', if for any < m < n and 
x > 0, there holds 

P< sup {\(n — m) — [a(n) — a(m)] } > x > < h(x). (10) 

0<m<n ' 

Through some manipulations, Eq.( fiTI]) can be expressed as the max-plus 
convolution: 

P{a(n) < a®\(n) - x} < h(x). (11) 

Here, a®A(n) can be considered as the expected time that the packet would 
arrive to the head-of-line (HOL) if the flow has passed through a virtual SSQ 
with the (deterministic) service curve A(n). The packet is expected to arrive 
not earlier than the expected HOL time. Here x represents the difference 
between the expected HOL time and the actual arrival time. The violation 
probability is bounded by the non-increasing function h(x). 

We use the v.w.d SAC to characterize the arrival traffic in Example [TJ 

Example 2. 

Consider a flow that consists of packets having the fixed packet size. 
Suppose that all packet inter-arrival times are exponentially distributed with 
mean -. Based on the steady-state probability mass function (PMF) of the 



queue-waiting time for an M/D/l queue [331 ]. we say that the flow has a v.w.d 
SAC \{n) = h- n with bounding function h exp for < h < -. Let p = u ■ h. 
We obtain the bounding function of the probability that the waiting delay 
W{n) exceeds x(> 0) 

LfJ 

h^(x) = 1 - (1 - p) J2 e -^ h - x ^ iH ~ x)]l 
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where, [y\ denotes the floor function. 

The definition of v.w.d SAC is more strict than that of i.a.t SAC. As a 
result, it is not trivial to derive the v.w.d SAC for an arrival process even 
if it can be characterized by an i.a.t SAC. Thus, it is important to explore 
whether there exists some relationship between the i.a.t SAC and the v.w.d 



Theorem 2. 1. If a flow has a v.w.d SAC X(n) G T with bounding func- 
tion h[x) G T , then the flow has an i.a.t SAC X(n) G T with the same 
bounding function h(x) G T . 
2. Conversely, if a flow has an i.a.t SAC X{n) G T with bounding function 
h(x) G Q, it also has a v.w.d SAC A_ ?? (n) G T with bounding function 



Remark. In the second part, h(x) G Q while not h(x) G T. If the 
requirement on the bounding function is relaxed to h(x) G the second 
part may not hold in general. 

Theorem H] reveals that if an arrival process can be modeled by a v.w.d 
SAC A(n), then X{n) is also the i.a.t SAC of this arrival process. On the 
other hand, if an arrival process can be modeled by an i.a.t SAC with the 
associated bounding function in Q, then this arrival process also has a v.w.d 
SAC which may be associated with a more loose bounding function. 

It is worth highlighting that the v.w.d SAC looks similar to the v.b.c SAC 
(see Definition [2]) defined in the space-domain. Since these two models play 
an important role in performance analysis in their respective domains, we 
establish their relationship in the following theorem. 

Theorem 3. 1. // a flow has a space-domain v.b.c SAC a(t) G J 7 with 
bounding function f(x) G T , the flow has a time-domain v.w.d SAC 
X(n) G T with bounding function h(y) G T , where 



SAC. 



h v (x) G Q , where for i] > 




J i 



A(n) = inf{r : a{r) > n}, and h{y) = 
with z _1 (y) denoting the inverse function ofy, where 
y = z{x) = sup{A(/c) — X{k — x)}. 

k>0 



2 Note that r\ should not be greater than lim. 



n 
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Specifically, if A(-) is sub- additive, z(x) = \(x). 
2. Conversely, if a flow has a time-domain v.w.d SAC X(n) G T with 
bounding function h(y) G T , the flow has a space-domain v.b.c SAC 
a(t) 6 J with bounding function f(x) G T , where 

a(t) = sup{fc : X(k) < t}, and f(x) = h^z^ 1 ^)) 

with z~ 1 {x) denoting the inverse function of x, where 

x = z(y) = sup{a(r + y) — a{r) + 1}. 

T>0 

Specifically, if a(-) is sub-additiv^, z(y) = a(y) + 1. 

Note that in Theorem [3j the arrival curve a(t) denotes the cumulative 
number of arrival packets while not the cumulative amount (in bits) of arrival 
traffic. 

The generalized stochastically bounded burstiness (gSBB) jssj] is a spe- 
cial case of the space-domain v.b.c SAC. A summarization of some well- 
known traffic belon ging to gSBB is given [22j, including both Gaussian self- 
similar processes ji[ U| 25] [30], such as fractional Brownian motion, and non- 
Gaussian self-similar processes, such as a— stable self-similar process [3] [24], 
and the (a(8), p(9)) stochastic traffic model 0)0]. With Theorem [31 the 
following example shows that gSBB can be readily represented using the 
time-domain v.w.d SAC. 

Example 3. 

If an arrival process A(t) can be described by gSBB with upper rate p 
and bounding function f(x) G i.e., for any t, x > 0, there holds 

P\ sup {A(s,t) - p ■ (t - s)} >x\ <f(x), 

then the process A(t) has a v.b.c SAC a{t) = p-t with the bounding function 
f(x). With Theorem [3] (1), the arrival process has a v.w.d SAC \(n) = - 
which is sub-additive and the bounding function h(y) = f(p-y), i.e., 



P 



{ sup {- • (n-m) - [a(n) - a(m)\\ > yX < f(p-y). 

^ 0<m<n P > 



3 [3] clarifies that a(t) defines a meaningful constraint only if it is subadditive. If a(t) 
is not subadditive, it can be replaced by its subadditive closure. 
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Remark. Theorem [3] allows us to readily utilize the results of gSBB 
traffic for time-domain models. If the traffic is more suitable for being char- 
acterized by the time-domain traffic models rather than the space-domain 
traffic models, then the transformation between two domains can facilitate 
the analysis. 

3.2. Time-domain Service Models 

Queueing theory characterizes the service process of a system based on 
the per customer service time. Like the arrival model, time-domain service 
models extend to the cumulative service time. 

If packet P(n) arrives to a network system after packet P(n — 1) has 
departed from the system, the departure time of P(n) is the arrival time 
a(n) plus the service time 5 n , i.e., a(n) + 5 n . If P(n) arrives to the system 
while P(n — 1) is still in the system, then its departure time is d(n — 1) + 5 n . 
The combination of both cases gives the departure time of P(n) 

d(n) = max[a(n), d(n — 1)] + S n (12) 

with d(0) = a(0) + 5q. Applying Eq. (JT2"j) iteratively to its right-hand side 
results in 

n 

d(n) = sup [a(m) + 5 k ] ■ (13) 

0<m<n , 
— — k=m 

The system usually allocates a minimum service rate to an arrival flow in 
order to meet its QoS requirements. The guaranteed minimum service rate 
is related to the guaranteed maximum service time for each packet of the flow. 
Accordingly, the time that the packet departs from the system is bounded. 
Denote the guaranteed maximum service time by 5 n . The Guaranteed Rate 



Clock (GRC) is defined based on S n [17| [18 



GRC(n) = max[o(n), GRC(n - 1)] + S n . (14) 

with GRC(0) = a(0) + 5 . Applying Eq. ffT4"|) iteratively to its right-hand side 
yields 

n 

GRC(n) = sup [a{m) + • ( 15 ) 

0<m<n , 
— — k=m 
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Eq.( TT5"j) is similar to Eq. (TT^|) except for that GRC{n) represents the guar- 
anteed departure tima3 while d(n) is the actual departure time. 

If Xir=m $i * s denoted by a function 7(71 — m + 1), then Eq. ffl5|) becomes 

GRC{n) = sup \a(m) + 7(n — m + 1)] = a®7(n) (16) 

0<m<n 



which is the basis for the time-domain (deterministic) service model [lOj. For 
systems that only provide service guarantees stochastically or applications 
that require only stochastic QoS guarantees, the service time may not need 
to be deterministically guaranteed. In this case, we extend the (deterministic) 
service curve into a probabilistic one. 

Definition 7. (i.d Stochastic Service Curve) 

A system is said to provide an inter-departure time (i.d) stochastic service 
curve 7(71) G T with bounding function j(x) G T , if for any n, x > 0, there 
holds 

pjd(n) - a<§)7(n) > x\ < j(x). (17) 

Note that the stochastic service curve of a service process is not unique. 
Therefore optimization is needed to find the SSC of a specific system. 

Example 4. 

Consider two nodes, the transmitter and the receiver. They communicate 
through an error-prone wireless link which is modeled as a slotted system. 
The wireless link can be considered as a stochastic server. Packets have 
fixed-length and are served in a FIFO manner by the transmitter. To sim- 
plify the analysis, we assume that the length of time slot equals one packet 



transmission timej 

The transmitter sends packets only at the beginning of a time slot. Due 
to the error-prone nature of the wireless link, the probability that a packet 
is successfully transmitted is determined by packet error rate (PER). Here, 
we assume that packet errors happen independently in every transmission 
with a fixed PER denoted by P e . The successful transmission probability of 



4 The guaranteed departure time is actually Gi?C(n)+error term |17j . where error term 
is determined by the employed service discipline. The underlying service discipline con- 
sidered throughout this paper is FIFO, under which, the error term is zero. 

5 It means we only compute the number of time slots in this example. 
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one packet is hence 1 — P e . If error happens, the unsuccessfully transmitted 
packet will be retransmitted in the next time slot immediately. In order 
to guarantee 100% reliability, the packet will be retransmitted until it is 
successfully received by the receiver. 

The per-packet service time 5 n is a geometric random variable with pa- 
rameter 1 — P e . The cumulative service time of successfully transmitting 
packets P(m) to P(n) is Xwc=m^ which follows the negative binomial dis- 
tribution with parameter 1 — P e . The mean service time denoted by 5 equals 



1-fV 

According to the complementary cumulative distribution function (CCDF) 
of the negative binomial distribution, the cumulative service time between 
two arbitrary packets P(m) and P(m + n) is given by 

p{ h > 6 ■ (n + 1) + x} < ( n J (1 ~ Pe) n+1 Pt (n+1) (18) 

for any x > 0, where |~-~| is the ceiling function. 

The right-hand side of Eq. ffl8l) represents the bound on the probability 
that the actual cumulative service time exceeds the cumulative mean service 
time. Let 7r,(ra) = 5 ■ n + r] ■ n for r\ > and j(x) denote the right-hand side 
of Eq. ffl8|) . From Definition [71 we know 

d{n) — a®7^(n) 

■n 

= sup [a(m) + Sk] — sup [a(m) + (5 + rf) ■ (n — m + 1)] 

0<m<n , 0<m<n 

— — k=m — — 



it 



— SU P [ ^fc — ^ " ( n — m + 1) — ' ( n — m + 1)] ) 

0<m<n , 
— — k=m 

from which, we have 

n 

P< sup [ 2, 5k — 5 ■ (n — m + 1) — 7] ■ (n — m + 1)] > x > 

L 0<m<n , ) 
— — k=m 

n n 

< y~] P| y~] 5 k -6 ■ (n-m + i) > x + ij- (n — m + l)j 

m=0 k=m 

n 

< \J j{x + r] ■ (n — m + 1)) 



m=0 
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n +l n oo 

^2j(x + r]- k) < - j(y)dy 
k=i 71 Jx 



Thus, we conclude that this error-prone wireless link provides an i.d SSC 
j v (n) with the bounding function j v (x) for r\ > 0, where 

j v (n) = 5 ■ n + r] ■ n and j v {x) = - I j{y)dy 



V Jx J i 

Since Eq. (1181) is only relevant to the cumulative service time and does not 
involve the arrival process, it provides a method to find the i.d SSC. 

Remark. Example H] demonstrates that we can obtain the i.d SSC from 
analyzing per-packet service time. However, if applying the space-domain 



results to this case, we need an impairment process 22] to characterize the 
cumulative amount of service consumed by unsuccessful transmissions. In 
other words, we still need to compute the cumulative slots due to failed 
transmission and then convert it into the amount of service. Such conversion 
may introduce error or result in looser bounds whereas the time-domain 
model directly computes the service time and avoids the conversion error. 
This simple example thus illustrates the feasibility of the time-domain service 
curve model. 

In Section HI we show that many results can be derived from the i.d 
SSC. However, without additional constraints, we have difficulty in proving 
the concatenation property for the i.d SSC. To address this difficulty, we 
introduce another service curve model in the following. 

Definition 8. (rj- Stochastic Service Curve) 

A system is said to provide an 77-stochastic service curve 7(72) G J 7 with 
bounding function j^x) G T , if for any n, x > 0, there holds 

P< sup \d(m) — a®7(m) — 77 • (n — m)] > x > < j v {x), (19) 

!> 0<m<n ' 

for any small rj > 0. 

Note that the left-hand side of Eq. (fT9l) represents a property that is typ- 
ically hard to calculate. It means that Definition [8] is more strict than Def- 
inition [71 Thus it is important to find the relationship between the i.d SSC 
and the //-stochastic service curves. 
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Theorem 4. 1. If a system provides to its arrival process an r\-stochastic 
service curve 7(71) with bounding function j v (x) G T , it provides to 
the arrival process an i.d SSC ^(n) with the same bounding function 
j n (x) G T; 

2. If a system provides to its arrival process an i.d SSC j(n) with bounding 
function j(x) G Q, it provides to the arrival process an rj-stochastic 
service curve ^(n) with bounding function j v {x) G Q for r] > 0, where 



3v\ x ) 



j{x) + - I j{y)dy 

r] 



Again, in the second part of Theorem HJ j{x) G Q while not j(x) G T . If 
the requirement on the bounding function is relaxed to j(x) G the above 
relationship may not hold in general. 

Definition [7] explores the relationship between the arrival process and the 
departure process, but it does not explicitly characterize the service process. 
From Eq. fri7|) . it is not trivial to find the stochastic service curve ^{n) for a 
specific system. Example H] illustrates how to add some increment r\ to the 
stochastic service curve 7(n). To this end, we expand Eq.(fl6l) to 

d{n) — a®^{n) = sup \a{m) + A(m, n)] — a®^(n). (20) 

0<m<n 

Without loss of generality, assume a(m ) (0 < m < n) is the beginning of 
the backlogged period in which packet P{n) is served. Then, 

sup \a{m) + A(m, n)\ = a(m ) + A(m , n) 

0<m<n 

and a®7(n) > a(mo) + j(n — m + 1). 

We rewrite the right-hand side of Eq. (120]) as 

a(m ) + A(m , n) — a®^(n) 
< a(m ) + A(m ,n) - a(m ) - 7(m + 1) 
= A(m Q ,n) — j(n — m + 1). (21) 

Note that Eq. fT2T]) holds for arbitrary m < n. Inspired by this, we define a 
new service curve model. 
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Definition 9. (Stochastic Strict Service Curve) 

A system is said to provide stochastic strict service curve 7(71) G T with 
bounding function j(x) G T ' , if the cumulative service time between two ar- 
bitrary packets P(m) and P(nf^ satisfies for any x > 0, 



Eq.( l2T|) reveals a relationship between the i.d SSC and the stochastic 
strict service curve. Furthermore, in Theorem H^2), the relationship between 
the stochastic strict service curve and the 77— stochastic service curve is ob- 
tained. 

Theorem 5. Consider a system providing stochastic strict service curve 
7(71) G J 7 with bounding function j(x) G T . 

1. It provides an i.d SSC 7(77) with the same bounding function j(x). 

2. If j(x) G Q , it provides an 77— stochastic service curve 7(77) with bound- 
ing function j v {x) G Q , where 



Note that the second part of Theorem [5] requires the bounding function 
j(x) G Q while not j(x) G T . 

4. Fundamental Properties 

In this section, we explore the four fundamental properties for time- 
domain models, i.e. service guarantees, output characterization, concate- 
nation property and superposition property. Some properties can only be 
proved for the combination of a specific traffic model and a specific service 
mode. This is why we have established various transformations between 
models in Section [3J With these transformations, we can flexibly apply the 
corresponding models to specific network scenarios. 

4-1- Service Guarantees 

Suppose that the arrival process has a v.w.d SAC and the service process 
has an i.d SSC. Under this condition, we derive the delay bound and backlog 
bound. 



6 If P(m) and P{n) are in the same backlogged period, A(m,n) = d(n) — d(m — 1). 




(22) 



j(x) + - 




J 1 
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4-1.1. Delay Bound 

The system delay significantly impacts QoS and is an important perfor- 
mance metric. 

Theorem 6. (System Delay Bound). 

Consider that a system provides an i.d SSC 7(72) G J 7 with bounding 
function j(x) £ J 7 to the input which has a v.w.d SAC X(n) G J 77 with 
bounding function h(x) G J 7 . Let D(n) = d(n) — a(n) be the system delay of 
packet P{n). For x > 0, D{n) is bounded by 

P{D(n) > x} < j <g> h([x - 7 A(l)] + ). (23) 



If the arrival process and the service process are independent of each 



other, we obtain another system delay bound according to Lemma 6.1 [22 



Lemma 1. (System delay bound: independent condition) 

Consider that a system provides an i.d SSC 7(n) G J 77 with bounding 
function j(x) G J 7 to the arrival process which has a v.w.d SAC A(n) G J 7 
with bounding function h(x) G T. Suppose that the arrival process and the 
service process are independent of each other. Then for x > 0, the system 
delay D{n) is bounded by 

P{D(n) > x} < 1 - ] * ^([a; - 7 A(l)] + ), (24) 

where j(x) = 1 — [j(a;)]i and h(x) = 1 — [/i(x)]i. 



4-1.2. Backlog Bound 

The system backlog represents the total number of packets in the system 
at time t, including both the packets waiting in the buffer and the packet 
being served. It is determined by function ([T]): 

B(t) < inf |/ > 0,sup{n > : a(n < t)} : d(n - I) < a(n)|. 

The following theorem provides a probabilistic bound on the system backlog 
for the given arrival process and service process. 

Theorem 7. (Backlog Bound) 

Consider that a system provides an i.d SSC j(n) G J 7 with bounding 
function j(x) G J 7 to the arrival process which has a v.w.d SAC X(n) G J 7 
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with bounding function h(x) £ J 7 . The system backlog at time t (> 0) is 
bounded by 

P{B{t) > x] < j ® h(^0X([x - 1]+)) (25) 



for x > 1 . 
Let 



H(X, 7 + x) = sup < inf [A; > : 7(771) + x < A(m + k)] > 

m>0 ^ > 



represent the maximum horizontal distance between functions X(n) and , y(n) + 
x. The probability that Bit) exceeds if (A, 7 + x) is bounded by 

p{B(t) > #(A,7 + x) + l} < j®h(x). (26) 

Remark. H(X,j+x) can be considered as the maximum system backlog 
in a (deterministic) virtual system, where the arrival process is A(n) and the 
service process is 7(n) + x. Eq. (l26|) is thus a bound on this maximum system 
backlog. 

If the arrival process and the service process are independent of each 
other, another backlog bound is derived according to Lemma 6.1 22 



Lemma 2. (Backlog Bound: independent condition) 

Consider that a system provides an i.d SSC ^(n) £ J 7 with bounding 
function j(x) £ J 7 to the arrival process which has a v.w.d SAC A(n) £ T 
with bounding function h(x) £ T . Suppose that the arrival process and the 
service process are independent of each other. Then the system backlog at 
time t (> 0) is bounded by: 

P{B(t) > x} < 1 - ] * h(-r®\([x - 1]+)) (27) 

for x > 1 . 

The probability that B{t) exceeds H(X, 7 + x) is bounded by 

P{B{t) > H(\,i + x) + 1} < l-j*h(x). (28) 

4-2. Output Characterization 

The previous section has presented how to derive the service guarantees 
in a single node. Another common scenario with which performance analy- 
sis deals is the end-to-end performance. An intuitive and simple approach 
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departure: d(n) 



Figure 1: Output characterization 



is called node-by-node analysis [19[ which requires characterization of the 
departure process from a single node. 

Let us consider a simple network as shown in Figure [TJ The departure 
process of Server 1 is the arrival process for Server 2. 

The delay bound in Server 1 can be derived from the result of Section 
14.1.11 To derive the delay bound in Server 2, we need to characterize the 
arrival process to Server 2, which is the departure process from Server 1. 
The problem is how to characterize the departure process from Server 1. 

Theorem 8. (Output Characterization) 

Consider that a system provides an i.d SSC 7(71) G J 7 with bounding 
function j(x) G T to its arrival process which has av.w.d SAC X(n) G T with 
bounding junction h(x) G T . The output has an i.a.t SAC A07(n — m — 1) 
with bounding function j®h{x) G T , i.e., for any < m < n — 1, there holds 

P |A07(n - m - 1) - [d(n) - d(m)] > x\ < j <S> h(x). (29) 

Remark. In Theorem[HJ the initial arrival process has a v.w.d SAC while 
the departure process has an i.a.t SAC. In order to derive the service guar- 
antees in Server 2, we need Theorem [2] (2) to transform the i.a.t SAC into a 
v.w.d SAC. Such transformation introduces a loose bounding function. The 
node-by-node analysis thus generates a loose end-to-end delay bound. Net- 
work calculus possesses an attractive property, concatenation property, which 
is used to deal with the end-to-end performance analysis. The comparison 
between the node-by-node analysis and the concatenation analysis reveals 
that the latter yields a tighter end-to-end delay bound [il| . 

The output characterization property however is very useful when ana- 
lyzing complicated network scenarios, such as Figure [2j where flows join or 
leave dynamically. In order to analyze the per-flow service guarantees, the 
departure process from each single node should be characterized using the 
arrival process to the node and the service process provided by the node. 
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Figure 2: Complicated network scenario 



Moreover, if the arrival process and the service process are independent 
of each other, the following lemma depicts the departure process. 

Lemma 3. ( Output Characterization: independent condition.) 

Consider that a system provides an i.d SSC 7(72) G J 7 with bounding 
function j(x) G T to its arrival process which has a v.w.d SAC A(n) G T 
with bounding junction h(x) G T . The output has an i.a.t SAC A* G T with 
bounding junction h* (x) G T , where 



4-3. Concatenation Property 

The concatenation property aims to use an equivalent system to represent 
a system of multiple servers connected in tandem if each server provides a 
service curve to its input. Then this equivalent system can be considered as 
a 'black box' which also provides the initial input with a service curve. 

In the following discussion, 7 fc and j k denote the stochastic service curve 
and bounding function of the kth server. For packet P{n), the time arriving 
to the kth server is a k (n) and the time departing from the kth server is d k {n). 
For a network of N tandem servers, the initial arrival is a(n) and the final 
departure is d{n). 

Theorem 9. (Concatenation Property) 

Consider a jlow passing through a system ojN nodes connected in tandem. 
If each node k(= 1,2,. ..,N) provides an i.d SSC 7 fc (n) G T with bounding 
function j k (x) G Q to its input, the system provides to the initial input a[n) 
an i.d SSC j(n) with bounding function j(x) , where 




(30) 



7(n 



7 1 ®7'® • • ■ ®7 ( Vi) 



n 




with 



7 ( fc fc _i),>) = 7 fc W + {k - 1) • V ■ n 
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for k = 2,...,N and rj > 0, and 

1 f°° 

j k >*(x) = \j h {x) + - J k (y)dy} 1 

Vk Jx 

for k = 1, .. . , N — 1 and rjk > . 

The proof of Theorem [9] utilizes the relationship between the i.d SSC and 
the //—stochastic service curve. The following lemma directly describes the 
service characterization of a system of nodes connected in tandem, where 
each single node provides an ??— stochastic service curve to its input. 

Lemma 4. Consider a flow passing through a system of N nodes connected 
in tandem. If each node k(= 1,2,.. .,N) provides an r\- stochastic service curve 
7 fc (n) G T with bounding function j k (x) G T to its input, i.e., 

P< sup {d k (m) - a k ®^ k (m) - rj ■ (n - m)} > x \ < j k {x), 

0<m<n ' 

then the system provides to the initial arrival process an i.d SSC j(n) with 
bounding function j(x): 

7 ( n ) = T 1 (§) 7 ^ . . . ®j N N-i)v( n ) 

where 7^_n^(w) = 7 fc (n) + {k — 1) • 77 • n, k = 2, N , for any small 77 > 0. 

Remark. The proof of the concatenation property reveals another reason 
of defining the 77— stochastic service curve model. 

4-4- Superposition Property 

The superposition property can be applied for multiplexing individual 
flows into an aggregated flow under the FIFO aggregate scheduling. The 
arrival process of the aggregate flow can be characterized by a stochastic 
arrival curve if the arrival process of each individual flow can be stochastically 
characterized by a stochastic arrival curve. Then we only need to analyze 
the service guarantees for the aggregate flow since all constituent flows are 
served equally. 
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4-4- 1- Superposition of Renewal Processes 

The superposition of multiple flows essentially falls into the research is- 
sue - superposition of renewal processes. In queueing networks, an individual 
server may receive inputs from different sources. It is reasonable to assume 
that the arrival process to a server is a superposition of statistically inde- 
pendent constituent processes [26j. The individual constituent processes are 
typically considered as renewal processes. A renewal process is a counting 
process in which the times between successive events are independent and 



identically distributed possibly with an arbitrary distribution [32 



The superposition of renewal processes has been widely studied since 



the original investigation by Cox and Smith [13j. However, the renewal 
property is not preserved under superposition except for Poisson sources. 
More precisely, the inter-arrival times in the superposition process become 
statistically dependent. This property cannot be captured by the renewal 



model 35 



In the following, we introduce how to characterize the superposition pro- 
cesses of multiple flows from a network calculus viewpoint. 

4 -4 -2. Arrival Time Determination 

First, we only consider the superposition of two flows denoted by F\ 
and i*2. Let «i(n), 02 (n) and a(n) be the arrival process of F±, F2 and the 
aggregate flow, respectively. As shown in Figure El F\ and F 2 are aggregated 
in the FIFO manner. If two or more than two packets which belong to 
different flows arrive simultaneously, they are inserted into the FIFO queue 
arbitrarily. 



Flow 1 
Flow 2 







-► 


FIFO 


Aggregate flow 
► 











Figure 3: Aggregation of two flows 

Figure HJdepicts that the arrival process of the aggregate flow is dependent 
on the arrival process of two constituent flows. 

Recall that P(n) denotes the (n + l)th packet of the aggregate flow. The 
same notation is also used for constituent flows F\ and F 2 . Thus, packet P(n) 
of the aggregate flow is either the mth packet of flow Fi (i.e., Pi{m — 1)) or 
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-► time 
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Figure 4: Packet arrival time 



the (n + 1 — m)th packet of flow i<2 (i.e., P2(ti — t«)), where < m < 71 + 1. 
When m = 0, it means no packet of flow i*\ arriving yet. When m = n + 1, it 
means no packet of flow F 2 arriving yet. By convention, we adopt a^n) = 
for n < 0. Since m takes value between to n+1, there are n+2 combinations. 

Theorem 10. Consider that two flows F\ and F 2 arrive to a network system 
and are aggregated into one flow Fa in the FIFO manner. Let ai(n) ; 02(71) 
and o(n) be the arrival process of flows F\, F 2 and F^, respectively. Then 
the packet arrival time of the aggregate flow is determined by 

a(n) = min < max [01(771 — 1), 02(71 — m)] > (31) 

0<m<n+l I L J J 

with 



o(0) = min j max [0, o 2 (0)] , max [ai(0), 0] | = min[ai(0), o 2 (0)]. 

We use an example to explain the underling concept of Theorem [lOj In 
Figure HI observe the arrival process of the aggregate flow at time t. Packet 
P(4) (arrival time: a(4) < t) is the last arrival packet, which is either packet 
Pi(m — 1) or packet P 2 {^ — 771), depending on which packet's arrival time is 
closer to time t, i.e., o(4) = max[oi(m — 1), a 2 (4 — m)] for m — 0, 1, 2, 3, 4, 5, 
the arrival time of packet P(4) is one element of the following set denoted 
by A 



A = |ai(4), a 2 (4), max[oi(0), o 2 (3)], max[ai(l), a 2 (2)], max[ai(2), o 2 (l)], 
max[ai(3), a 2 (0)] \, 
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i.e., a(4) = min{A}. We notice that min{A} is actually the expansion of 
Eg. (13 II) . According to the packet arrival times of two constituent flows shown 
in Figure HI we have 

a(4) = min |ai(4), a 2 (4), max[ai(0), a 2 (3)] = a 2 (3), max[ai(l), a 2 (2)] = a 2 (2), 
max[ai(2), a 2 (l)] = ai(2), max[ai(3), a 2 (0)] = Oi(3)| = ai(2), 

which is consistent with Figure HI 

Theorem [TU1 can be generalized to the aggregation of N(> 2) flows. 

Corollary 1. Consider that N(> 2) flows F\,F2,...,Fn arrive to a network 
system and are aggregated into one flow Fa in the FIFO manner. Let ai(n), 
a 2 (n) ; ... ; ajv(n) and a(n) be the arrival process of the N constituent flows and 
the aggregate flow, respectively. Then the packet arrival time of the aggregate 
flow is determined by 

a(n) = min < max[ai(mi — 1), a 2 (m 2 — 1), ajv(ra 

y] mi=n+l,mjS:\0,n+l\ L 

with 

a(0) = min {ai(0),a 2 (0),...,ajv(0)}. 

4-4-3 ■ Superposition Process Characterization 

Eq. (l3~Tj) can compute the packet arrival time of the aggregate flow. How- 
ever, we still have the difficulty in characterizing the packet inter-arrival time 
of the aggregate flow if the packet inter-arrival times of two constituent flows 
follow the general distribution. For this reason, it is difficult to directly 
characterize the arrival process of the aggregate flow from the temporal per- 
spective. Alternatively, we rely on the available results of the superposition 
property explored in the space-domain (see Theorem [TJ. 

In the space-domain, the traffic arrival process is characterized based on 
the cumulative amount of arrival traffic. In the following, we use A(t), Ai(t) 
and »4. 2 (t) to denote the cumulative number of arrival packets of the aggregate 
flow up to time t, the cumulative number of arrival packets of F\ up to time t 
and the cumulative number of arrival packets of F 2 up to time t, respectively. 
A{t) is the sum of A\{t) and *4 2 (£), from which we can find the stochastic 
arrival curve for the aggregate flow. 



N-l 
i=l 

(32) 
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Condition Target 




Figure 5: Transformation in Theorem 1111 

As shown in Figure El the condition is that the time-domain stochastic 
arrival curve of all constituent flows are known, and the target is to verify 
that the aggregate flow also has a time-domain stochastic arrival curve. 

If a flow has a time-domain v.w.d SAC, with Theorem [3(2), this flow has 
a space-domain v.b.c SAC, for which the superposition property holds (refer 
Theorem [T]). Applying Theorem [3(1) gives rise to the v.w.d SAC for the 
aggregate flow. 

If flow Fi has a v.w.d SAC Aj(n) with bounding function hi(x), i = 
1,2,...,N, from Theorem [3(2), we can verify that flow F$ has a v.b.c SAC 
ai(t) with bounding function fi(x) = hAzT (x)\ where cti(t) and z7 (x) are 
given in Theorem [3(2). Furthermore, according to Theorem [TJ, the aggregate 
flow has a v.b.c SAC ait) = with bounding function f{x) = 

fi g) • • • (g) /jv(x). Finally, we apply Theorem [3(1) and can verify that the 
aggregate flow also has a v.w.d SAC. 

Theorem 11. (Superposition property) 

Consider the aggregate of N flows. If the arrival process of each flow has 
a v.w.d SAC Xi(n) E F for i = 1, 2, N, i.e., 

P{ai{n) < ai®\i(ri) - y} < h^y), 

which implies that every flow also has a v.b.c SAC 

ai{t) = sup{/c : \i{k) < t} 

with bounding function 

f i (x) = h i {zr\x)) 
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where z i (x) denote the inverse function of x: 

x = Zi(y) = sup{a,(r + y)- a^r) + 1}. 

r>0 

Then the aggregate arrival process a(n) has a v.w.d SAC \(n) with bounding 
function h(y), where 

N 

A(n) = inf{r : £>(r) > n}, h(y) = /(^(y)), 
i=i 

/(x) = /i £g> ■ • ■ <S> /at(^) awd denoting the inverse function of y: 

y = z(x) = sup{A(A;) — X(k — x)}. 

k>0 

4-4-4- Special Case: Superposition of Poisson Processes 

As we have mentioned in Section 14.4. H the Poisson process is a special 
case of renewal processes because its renewal property is preserved under 
superposition. In addition, the superposition of multiple Poisson processes is 
still a Poisson process. From the temporal perspective, the inter-arrival time 
between two arbitrary events of a superposition of Poisson arrivals follows 
the Gamma distribution. 

Example 5. Consider the superposition process of two independent Poisson 
arrival processes. Suppose that all packets of both arrival processes have the 
same size. The packet inter-arrival times of two Poisson processes follow 
exponential distributions with mean — j- and respectively. Find the time- 
domain v.w.d SAC for the superposition process. 

In Example [2] the v.w.d stochastic arrival curve for a Gamma process has 
been derived. We thus know that the superposition process has a v.w.d SAC 
X s (n) = T s ■ n (0 < T s < ) with bounding function h s (x): 

I— J 

K{x) = i - (i - Ps ) jr e -&»+*x.T.-«) [(ai^Zkzf)!!, 

i=0 

where p s = (/ii + fi 2 ) ■ T s . 

Remark. It is readily to generalize the above example into the superpo- 
sition of multiple independent Poisson processes. 
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5. Conclusions and Open Issue 



This paper presented a temporal network calculus to formulate queueing 
systems in communication networks where applications can tolerate a certain 
level of performance violation. The time-domain models make it feasible 
to characterize the temporal behavior of network traffic and capture the 
temporal nature of the network capacity perceived by individual packets. 

The models are defined in such a way to compromise between simple 
models and complex models. The former may not be sufficient to explore the 
fundamental properties whereas the latter may be too difficult to build. In 
order to solve this dilemma, we propose a transformation method such that 
the appropriate models are selected to some specific scenario. Moreover, 
we also link the temporal network calculus and the existing space-domain 
network calculus results through connecting the time-domain v.w.d arrival 
curve and the corresponding space-domain v.b.c arrival curve. 

Four properties investigated in the time-domain facilitate performance 
analysis of various network scenarios. In addition, the proof of the super- 
position property has given insights into the importance of model transfor- 
mation. We believe that this temporal network calculus is applicable for 
analyzing networks where users are served probabilistically and compliments 
the current network calculus results. 

The leftover service characterization is useful for per-flow performance 
analysis and has been proved in the space-domain. We attempted to tackle 
this property under the condition that the arrival process has a deterministic 



time-domain arrival curve and the service process provides an i.d SSC [36 
We will expand the investigation to this property under the general condition. 
One challenge is to be able to decouple the constituent flow's arrival process 
from the aggregate arrival process. 
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Appendix A. Proofs of theorems and lemmas 

Proof of Theorem [2] 

Proof. The first part follows from that for any < m < n, there trivially 
holds 

\(n — m) — \a(n) — a{m)] < sup {A(t7 — m) — \a{n) — 0,(777)]}. 

0<m<n 

For the second part, there holds 



P< sup {A_ ?? (77 — m) — [0(77) — a(m)]} > x 

0<m<n 

< P< sup {A_,(n - m) - [a(n) - fl(m)]} + > 1 . 

^ 0<m<n ' 

For any x > 0, 

p|{A(t7 — m) — 77 • (n — m) — [a(n) — a(m)]} + > x 

= P |a(t7 — m) — 77 • (77 — 777) — [0.(77) — a(m)] > x 

= p|a(?7 — m) — [a{n) — a(m)] > a; + 77 ■ (n — m) 
< h[x + 77 • (77 — 777)) . 

Based on the above steps, we have 



Pi sup {A_^(t7 — m) — [a{n) — a(m)]} > x 

L 0<m<n 

n 

< p|{A_^(77 - 777) - [0.(77.) - a(777)]} + > x j 

m=0 
n 

< + 77 • (77 — m)) 

m=0 
n 

= /i(a; + i] ■ k) 

k=0 

00 

< + i] ■ k) 



k=0 
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k=l 



< /i(a?) + - I h{y)dy. 

V Jx 

The right-hand side of the last inequality still belongs to Q. The second part 
follows from the above inequality and the fact that the probability is always 
not greater than 1. □ 

Proof of Theorem [3] 

Proof. (1) From Lemma 2 [36], we know that for any t,x > 0, event 

{A(t) < A® a(t) +x} 

implies event 

{a(n) > a<8)\(n) — y} 

where 



y = sup {X(k) — X(k — x)} = z(x). 

k>0 



Thus, there holds 



P{A(t) < A®a(t) + x} < P{a(n)>a®X(n)-y} 
=^> P{a{n) < a ®\(n)-y) < P{A{t)> A®a(t)+x). 

< f(x), 

Particularly, if A is sub-additive, i.e. X(a + b) < X(a) + X(b) for any a and b, 
we then have: 

P{a(n) < a®X(n) - X(x)} 

< P{a(n) < a®X(n) - sup[A(A;) - X(k - x)] } 

k>0 

< f{x). 

Hence, the first part follows. 

(2) From Lemma 3 36], we know that for any n, y > 0, event 

{a(n) > a®X(n) — y} 
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implies event 

{A(t) <A® a(t) + x} 

where 

x = sup{a(u + y) — a(u) + 1} = z(y). 

u>0 

Thus, there holds 

P{a(n) > a®\(n) - y} < P{A(t) < A <g> a{t) + x) 

=> P{A(t) > A® a(t) + x) < P{a{n) <a®X{n) -y) 

< Ky). 

Particularly, if a is sub-additive, we have 

P{A(t) > A ® a(t) + a(y) + 1} 

< P{A(t) > A ® a(t) + sup[a(u + y) - a(u) + 1] } 

< %), 

which ends the proof. 
Proof of Theorem [4] 

Proof. The first part follows since there always holds 

d{n) — a®7(n) < sup [d(m) — a®7(m) — r] • (n — m)l 

0<m<n 

by letting m = n on the right hand side. 
For the second part, there holds 

sup \d{m) — a®7(m) — r) • (n — m)] 

0<m<n 

< sup {d(m) — a®7(m) — 77 • (n — m)} + . 

0<m<n 

Hence for any x > 0, there exists 

P{ sup — a<g>j(m) — 77 • (n — m)} > 2} 

0<m<n 

n 

< P\d{m) — a®7(m) — r] • (n — m) > x} 

m=0 

36 



n 
u=0 

< \j{x) + - j{y)dy 

L VJx 

The right-hand side of the above inequality still belongs to Q and is always 
not greater than 1. The proof of the second part is completed. □ 

Proof of Theorem [6] 

Proof. For any n > 0, according to the definition of D(n), there holds 

D(n) = d(ri) — a(n) 

= [d(n) — a®7(?i)] + [a®7(n) — a(n)] 

= [d(n) — a®7(ri)l + sup |a(m) + j(n — m + 1) — a{n) j 

0<m<ra 

= \d{n) — 0,07(72)] + sup {A(n — m) — \a{n) — a{m)} + 



0<m<n 



7(71 — m + 1) — A(n — m)} 

< <i(n) — a®7(n) + sup {\(n — m) — \a{n) — a(m)]} 

0<m<n 

+ sup {7(71 — m + 1) — X(n — m)} 

0<m<n 

< d(n) — a(g)7(n) + sup {A(t7 — m) — [0,(77) — a(m)]} 

0<m<n 

+ sup{7(fc + 1) - A(fc)}. 

fc>0 

To ensure system stability, we require 

lim -r[j(k) - \(k)] < 0. (A.l) 

In the proofs of the following theorems, without explicitly stating, we shall 
assume Eq. (lA.lj) holds. 

In addition, the following results are given 

P{d(n) — 007(77) > a;} and P< sup {A(n — m) — [a(n) — o(m)] } > x >. 

^ 0<m<n ' > 

From Lemma 1.5 [22J and sup fc>0 {^(k + 1) — A(A;)} =70 A(l), we can 
conclude 

P{D(n) > x} < j ®h(x - 7 A(l)). 
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□ 

Proof of Theorem [7] 

Proof. According to the backlog definition 

B{t) < inf jx > 0,sup{n > : a(n) < t} : d(n - x) < o(n)}, 

we need to prove the bounding function on the violation probability, i.e., 
P{B(t) > x}. For ease of exposition, let n = m + x, then we have 

d{m) — aim + x) 
= \d(m) — 0(8)7 (m)l + \a®^im) — a{m + x)l 
= \d(m) — a(§)7(m)l + sup {a(k) + 7(777. — k + 1)} — aim + x) 

0<fc<m 

= [d(m) — a(8>7(m)] + sup |A(m + x — /c) — [aim + x) — a(/c)] 

0<fc<m 

+ 7(771 — + 1) — A(m + x — fc)} 

< [d(m) — 0(8)7(777.) 1 + sup {A(m + x — fc) — [a(m + x) — a(fc)]j 

0<k<m+x 

— inf |A(m — + x) — 7(777. — fc + 1)} 

0<k<m 

Let v = m — k + 1. The above inequality is written as 
d(m) — a(m + x) 

< \d(m) — a®j(m)] + sup {A(m + x — fc) — [a(m + x) — a(fc)]j 

0<fc<m+i£ 

~i^ f ^ {A(w + x-l) - 7 (w)}. 

l<^<m+l 

Because there holds 

inf {A(f + x — 1) — 7(f)} > inf |A(tj + x — 1) — 7(f)} 

l<v<m+l v>l 

= A0 7 ([x-1] + ), 

with the same conditions as analyzing the delay, we obtain 
P{B(t) >x}<j® /7(A0 7 ([x - !]+)). 
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To prove Eg. (125]) . we replace x = H(X,j + y) + 1 in event {B{t) > x} and 
have 

d(m) - a(m + H(X, 7 + y) + 1) 
< [d(m) - a®j(m)] + a(§X (m + H(X, 7 + y) + l) 

-a(m + #(A, 7 + y) + l) + sup {7^) - \(v + H(X, 7 + y))}. 

v>0 

The definition of H(X, 7 + 7/) implies 

7(«)+i/< A(« + #(A,7 + jO) 

for any t> > 0, i.e., 

sup {7O) - X(v + H(X, 7 + y)) } <-y. 

v>0 

Then we conclude 



□ 



P{B{t) > H(X, 7 + x) + 1} < j ® /i(x). 

Proof of Theorem [8] 

Proof. For any two departure packets m < n, there holds 

d(m) — <i(n) < <i(m) — a(n) 

= d{m) — a{n) + a®^{m) — a®7(m) 

= [c?(m) — a®7(m)] + sup {a(k) + 7(m — k + 1)} — a(rz) 

0<fc<m 

= [d(m) — a®7(m)] + sup {7(771 — k + 1) — [a(n) — a.(A:)]} 

0<fc<m 

= \d(m) — a®7(m)] + sup < 7(777 — k + 1) — A(77 — /c) 

+ A(t7 — k) — [a(n) — a(k)] 
< \d(m) — 0(817(777)] + sup < A(r7 — k) — [a(n) — a(k)] 

0<k<m < 

+ sup {7(771 — k + 1) — A(ti — k)}. 

0<k<m 
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Let v — m — k + 1. Then the above inequality is written as 

d(m) — d(n) < \d(m) — 0(87(777)] + sup < X(n — k) — \a(n) — a(k)] > 

o /••• i, J 

— inf {A(t7 — m — 1 + v ) — 7(f)} 

l<v<m+l 

< [d(m) — 0(87(777)] + sup < A(n — k) — \a(n) — a(k)} > 

0<k<n ^ > 



- inf { A(n - m - 1 + - 7 (u) } 

where the last step is because 

inf [/J < inf [/J. 

Adding inf <„< m +i {A(n — 777 — 1 + v) — 7(f)} to both sides of the above 
inequality results in 

inf {A(n — m — 1 + v) — 7 (f)j — [d(n) — d{m)] 

0<v<m+l 

< \d{m) — 0(8)7(771)] + sup {A(n — fc) — [0(77) — a(fc)] }. 

0<fe<n 

In addition, there holds 

A0 7 (n — m — 1) = inf {A(n — m — 1 + i> ) — 7 (t> )} 

< inf {A(n — m — 1 + v ) — 7 (v)j. 

0<t)<m+l 

To ensure that the right-hand side of the above inequality is meaningful, it 
requires n — m — 1 > 0. With the same conditions as analyzing delay, we 
conclude 

p|A0 7 (n - 777 - 1) - [d(n) - d(m)\ > x| 

< P< [d(m) — 0(87(777)] + sup {\(n — k) — [a(n) — a(k)]} > x > 

0<k<n ' 

< j<g>h(x). 

□ 
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Proof of Theorem [fJ] 

Proof. We shall only prove the three-node case, from which, the proof can 
be easily extended to the iV-node case. The departure of the first node is 
the arrival to the second node, so d}{n) = a 2 (n) and d 2 (n) = a 3 (n). We then 
have, 

d{n) - a®7 1 (t)7j(8)72 J? (n) 

= d(n) — sup < a®7 1 (m) + 7^®72„(n — m + 1) [ + <i 1 (m) — d^m) 
0<m<n l> '' J 

< d(n) — sup < 7^®7L( n — to + 1) + d 1 {m) — 77 ■ (n — to + 1) 

0<m<n 

— \d l {m) — a®7 1 (m) — rj • (n — to)] j 

< d(n) — sup < 7^®72 T? («- — m + 1) + d l (m) — 77 ■ (n — to + 1) f 

0<m<n -* 

+ sup < d 1 (m) — a®7 1 (m) — 77 ■ (n — m)] > 

0<m<n ^ * 

= d(n) — sup {a 2 (m) + sup [7 2 (A;) + 77 • A; + 7 3 (77, — m + 1 — A;) 

0<m<n 0<fc<n— m+1 

+2r] ■ (n — m + 1 — k)) — rj ■ (n — m + 1) } + 
+ sup {^(m) — a®7 1 (m) — 77 • (n — to)} 

0<m<n 

= d(n) — sup |a 2 (m) + sup [7 2 (A;) + 7 3 (n — to + 1 — A;) 

0<m<n 0<fc<n— m+1 

+77 • {n — to + 1 — A;)]} + sup {d(m) — a®7 1 (m) — 77 • (n — to)} 

0<m<n 

= d{n) — a 2 ®7 2 £S>7 3 (n) + sup |c? 1 (m) — a®7 1 (m) — 77 • (n — to)} 

0<m<n 

< c?(n) — sup < a 2 ®7 2 (m) + 7 3 (tt, — to + 1) f — a 3 (m) + 77 ■ (n — to + 1) 
+d 2 (m) — 77 • (n — to) + sup {(^(m) — 0(8)7 (to) — 77 • (n — to)} 

0<m<n 

< d{n) — sup |a 3 (TO) + 7 3 (n — to + 1) — 77 • (n — to + 1)} 

0<m<n 

+ sup |c/ 2 (to) — a 2 ®7 2 (m,) — 77 ■ (n — to)} 

0<m<n 

+ sup {^(to) — o®7 1 (to) — 77 • (n — to)} 

0<m<n 
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d{n) — a 3 (8>7 3 (n) + sup {d 2 (m) — a 2 cg)7 2 (m) — rj ■ (n — m)} 

0<m<n 

+ sup |<i 1 (m) — a®7 1 (m) — i] ■ [n — m)}. 



0<m<n 

Based on the relationship between the i.d SSC and the ^-stochastic service 
curve presented in Theorem H(2), the following inequality holds 

P{d(n) - a®7 1 ®7 2 ®7^(n) > x} < f <g> f m ® j 1 ' 111 , 

which completes the proof. 

Note that both the max-plus convolution and the min-plus convolution 
are associative and commutative. □ 

Proof of Lemma 3] 

Proof. We shall only prove two-node case, from which, the proof can be 
extended to the iV-node case. Keep in mind that a 2 (n) = d 1 (n). For the 
two-node case, we have 

d(n) - a®^ 1 ^ 2 ^) 
= d{n) — sup |a(g)7 1 (m) + 7 2 (n — m + 1) + 77 • {n — m + 1)} 

0<m<n 

< d(n) — sup |a®7 1 (m) + 7 2 (n — m + 1) + 77 • [n — m)} 

0<m<n 

+d 1 (m) — a 2 (m) 

< d{n) — sup |a 2 (m) + -y 2 (n — m + 1)} 

0<m<n 

+ sup |<i 1 (m) — a®^ 1 (rri) — ?7 • {n — m)} 

0<m<n 

= d{n) — a 2 ®7 2 (n) + sup |<i 1 (m) — a®7 1 (m) — 77 ■ {n — m)} 

0<m<n 

< sup |<i(m) — a 2 Cg>7 2 (m) — 77 • (77, — m)} 

0<m<n 

+ sup {d 1 (m) — a®7 1 (^i) — ?7 • (n — m)}. 

0<m<n 

The last step holds because of Theorem H](l). From the condition, we con- 
clude 

P{d(n) - a®7 1 «)7j(n) > x] < j 1 <g> j 2 (x). 

□ 
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Proof of Theorem [TU1 

Proof. We use the induction way to prove this theorem. 

Step (1) We start from n = 1 with the given condition a(0) = min[ai(0), 02(0)]. 
If a(0) = cti(O), then a(l) = min jai(l), 02(0)}; if a(0) = 012(0), then 
a(l) = min |ai(0), a 2 (l)}. 

We expand Eq. fl3Tj) into the following expression: 



a(l) = min{ai(l),a 2 (l),max[ai(0),a 2 (0)]} 

un {«a(l), a 2 (0)}, if min[ai(0), a 2 (0)] = ai(0); 
min |ai(0), a 2 (l)}, if min[ai(0), a 2 (0)] = a 2 (0). 



Thus Eq. flBTj) holds for n — 1. 

Step (2) Assume n = k holds for > 1: 

a(fc) = min < max [ai(m — l),a 2 (/c — m)l \ (induction hypothesis), 

0<m<fc+l L J 

which has four solutions as below: 



a(k) 



f ai(k), if ai(k) < a 2 (0); 

a 2 (k), if a 2 (k) < ai(0); 

ai(m* — 1), if a 2 (/c — m*) < a\{m* — 1) for < m* < k + 1; 

^a 2 (A;-m*), if oi(m* — 1) < a 2 (A; — m*) for < m* < k + 1. 



Step (3) Prove n = k + 1 holds: 

a(fc + 1) = min < max [ai(m — 1), a 2 (k + 1 — m)l ^ 

0<m<(fe+l)+l I L J J 

which can be expanded into 

a(k + 1) = min |ai(A; + 1), a 2 {k + 1), max [ai(0), a 2 (&)] , max [cti(l), a 2 {k — 1)] , 
max [ai(2), a 2 {k — 2)] , max jai(A;), a 2 (0)] |. (A.2) 

We prove Eq. (IA.2p based on the four solutions of the induction hypothesis, 
respectively. 
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(I) If a(k) = ai(k) which implies a>i(k) < 02(0), then a(k + 1) is either 
a x (A; + 1) or a 2 (0) depending on the minimum one, i.e., a{k + 1) = 
min {ai(k + l),a 2 (0)}. Since the condition a\(k) < a 2 (0) implies 

ai(0) < ai(l) < ... < a x {k) < a 2 (0) < a 2 (l)... < a 2 (k), 
the induction hypothesis given in Step (2) is expanded into 

a(k) = min {a,i(k), a 2 (k), a 2 (k — 1), a 2 (0)} = ai{k). 

From this, Eq. (IA.2j) becomes 

a(k + 1) = min{ai(A; + 1), a 2 (k + 1), a 2 (k), a 2 (k — 1), a 2 (0)} 
= min {ai(/c + 1), a 2 (0)}, 

which proves that Eq. (lA.2p holds. 
(II) If a(A;) = a 2 (k) which implies a 2 (k) < ai(0), then a(k + 1) is either 
a 2 (k + 1) or ai(0) depending on the minimum one, i.e., a(k + 1) = 
min [a 2 (k + l),ai(0)}. Since the condition a 2 (k) < ai(0) implies 

a 2 (0) < a 2 (l) < ... < a 2 (k) < a^O) < ai(l)... < ai(k), 
the induction hypothesis given in Step (2) is expanded into 

a{k) = min {a 2 (k), ai(fc), a\(k — 1), ai(0)} = a 2 {k). 

From this, Eq. f ]A.2j) becomes 

a{k + 1) = min |a 2 (A; + 1), ai(/c + 1), ai(fc), ai(fc — 1), ai(0) j 
= min {a 2 (/c + 1), ai(0)}, 

which proves that Eq. flA.2l) holds. 
(Ill) Without loss of generality, if a(k) = ai(m* — 1) for < m* < k + 1, 
which implies a 2 (k — m*) < ai(m* — 1), then a(k + 1) is either ai(m*) 
or a 2 (k — m* + 1) depending on the minimum one, i.e., a(k + 1) = 
min |ai(m*), a 2 (A;+l— m*)}. Since the condition a 2 (k—m*) < ai(m* — 
1) implies a 2 (k — m*) < ai(m* — 1) < ai(m*) < a\{k) and a 2 (k — m*) < 
ai(m* — 1) < a 2 (k + 1 — m*) < a 2 (A;) for < m* < A;+ 1, the induction 
hypothesis given in Step (2) is expanded into 

a(k) = min \ai(k), a 2 (k), a 2 (k — 1), a 2 (k + 1 — m*), ai(m* — 1), 
ai(fc — 1)} = ax(m* — 1) 
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From this, Eq.( lA.2p becomes 



a(k + 1) = min {ai(k + 1), a 2 (k + 1), a 2 (k), a 2 {k — 1), 

a 2 {k + 1 — m*), a^m*), ai(m* + 1), ai(fc)} 
= min {a 2 (k + 1 — m*), ai(m*)} 

which proves that Eg. (1A.2I) holds. 
(IV) Without loss of generality, if a(k) = a 2 (k — m*) for < m* < k + 
1, which implies ai(m* — 1) < a 2 (k — m*), then a(k + 1) is either 
ai(m*) or a2(/c — m* + 1) depending on the minimum one, i.e., a(k + 
1) = min {ai(m*), a 2 {k + 1 — m*)}. Since the condition ai(m* — 1) < 
a 2 (k — m*) implies ai(m* — 1) < a 2 [k — m*) < a!(m*) < ai(k) and 
ai(m* - 1) < a 2 (k-m*) < a 2 (k + l-m*) < a 2 (k) for < m* < k + 1, 
the induction hypothesis given in Step (2) is expanded into 

a{k) = mm {ai(k),a 2 (k),a 2 (k — 1), ...,a 2 (k — m*),ai(m*), 



which proves that Eq. flA.21) holds. 
Combining the above three steps concludes that Eq. (1311) holds for all n > 



ai(k — 1)} = a 2 {k 



— m 



From this, Eq. (1A.2I) becomes 



a(k + 1) 



min [ai(k + 1), a 2 (k + 1), a 2 (k), a 2 (k - 1), 
a 2 {k + 1 — m*), ai(m*), ai(m* + 1), ai(fc) j 
min {a 2 (k + 1 — m*), ai(m*)} 



0. 



□ 
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